Distributions of Conductance and Shot Noise and Associated Phase Transitions 



OO 
O 

o 

CD 

in 



CD 



I 

O 

o 



> 

0^ 
On 
p 

a^ 

O 
OO 
O 



X 



Picrpaolo Vivo^, Satya N. Majumdar^ and Oriol Bohigas^ 
^ Department of Mathematical Sciences, Brunei University, 
Uxbridge, Middlesex, UBS 3PH (United Kingdom) 
^ Laboratoire de Physique Theorique et Modeles Statistiques (UMR 8626 du CNRS), 
Umversite Paris-Sud, Bdtiment 100, 91405 Orsay Cedex (France) 
(Dated: September 5, 2008) 

For a chaotic cavity with two indentical leads each supporting channels, we compute ana- 
lytically, for large A'', the full distribution of the conductance and the shot noise power and show 
that in both cases there is a central Gaussian region flanked on both sides by non-Gaussian tails. 
The distribution is weakly singular at the junction of Gaussian and non-Gaussian regimes, a direct 
consequence of two phase transitions in an associated Coulomb gas problem. 
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Ballistic transport of electrons through a mesoscopic 
cavity such as a quantum dot has been studied intensively 
both theoretically and experimentally in recent times [1] . 
In the simplest setting one considers a single cavity of 
submicron dimension (e.g., a quantum dot) with two 
identical leads (each supporting N channels) connecting 
it to two separate electron reservoirs. An electron, in- 
jected through one lead, gets scattered in the cavity and 
leaves by either of the two leads. The transport of elec- 
trons through such an open quantum system is encoded 

fr t' 

in the 2N x 2N unitary scattering matrix "5* = ( ^ / 

connecting the incoming and outgoing electron wavefunc- 
tions, where r, t are N x N reflection and transmission 
matrices from the left and r',t' from the right. Several 
experimentally measurable transport observables such as 
the conductance (i.e., the time averaged current) and the 
shot noise power (that describes the current fluctuations 
associated with the granularity of electronic charge e) can 
be expressed in terms of N transmission eigenvalues Ti 's 
of the N X N matrix T = tt'' . For example, the dimen- 
sionless conductance G ~ Tr(ttt) — X^iLi ^« [-] ^'^'^ ^^le 
shot noise P = Ti[tt^l~tt^)] = J^Zi T^{l~Ti) [3]. The 
eigenvalue < T,; < 1 has a simple physical interpreta- 
tion as the probability that an electron gets transmitted 
through the z-th channel. 

Over the past two decades, the random matrix theory 
(RMT) has been successfully used [1, 4] to model the 
transport through such a cavity. Within this RMT ap- 
proach, one assigns a uniform probability density to all 
S matrices belonging to the unitary group: this amounts 
to drawing S at random from one of Dyson's Circular 
Ensembles [5]. This, in turn, induces a probability mea- 
sure over the transmission eigenvalues Tj's whose joint 
probability density (jpd) is known [1, 4, 6] 
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constant and the Dyson index (3 characterizes different 
symmetry classes (/3 = 1, 2 according to the presence or 
absence of time-reversal symmetry and (3 = A 'm case of 
spin-flip symmetry). Given this jpd in (1) one is then 
naturally interested in computing the statistics of trans- 
port observables such as the conductance G and the shot 
noise P. 

The average and variance of G and P are known [1, 3] 
and in particular for large iV > 1, (G) = N/2, (P) = 
iV/8, var(G) = 1/8/3, var(P) = 1/64/3. The fact that the 
variance of G becomes independent of N for large N has 
been dubbed the 'universal conductance fluctuations' [1]. 
In contrast, much less is known for their full distributions. 
For the conductance, the full distribution is known ex- 
plicitly only for A'' = 1 and iV = 2 [4] and very recently, a 
solution in terms of Painleve trascendents for the special 
case /? = 2 has been announced [7]. Also, the behavior of 
these distributions very close to the endpoints (e.g., near 
G = and G = N) have been computed recently [8]. 

The purpose of this Letter is to present exact results 
for the fuU distributions V{G,N) and r{P,N) for all /3 
in the large A^ limit. For the conductance G, we show 

that for large A^, V{G,N) « cxp (^-fA^^^c (§)) where 
w stands for the precise asymptotic law: 
2\tlV{Nx,N)' 
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where = /o^ ■ • P {{Ti}) Y\, dTi is a normalization lar large A^ behavior: V{P, N) w exp [~%N^^p 



and the rate function g{x) is explicitly given in (10). 
Thus the distribution V{G, N) has a pure Gaussian form 
around the mean (G) — N/2 over the region A^/4 < G < 
3A^/4 and outside this central zone, it has non-Gaussian 
large deviation tails. The distribution has an extraor- 
dinarily weak singularity at the points G = A/4 and 
G = 2>N/A (only the 3-rd derivative is discontinuous!). 
One of our central results is to show that these two weak 
singularities in the conductance distribution are the di- 
rect consequence of two phase transitions in an associated 
Coulomb gas problem. The distribution of the shot noise 
P in the allowed range < P < A/4 displays a simi- 
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where \E'p(x) is given in (12). Once again the cen- 
tral Gaussian regime around the mean (P) = iV/8 is 
flanked by two non-Gaussian tails with a weak third or- 
der singularities at the transition points P = N/16 and 
P = 3A^/16. In addition to these distributions, we have 
also computed exactly the large N statistics of the inte- 
ger moments %i = X^t^i "^i"- Even though our analytical 
results are valid for large N, we find, rather remarkably, 
that the numerical results for relatively small values of 
N (such as = 4) agree fairly well with our asymptotic 
results. 

We start by demonstrating how the computation of 
'P{G,N) can be mapped to the calculation of the par- 
tition function of a Coulomb gas problem. By def- 
inition, V{G,N) = (<5(G-^^T,)) where () denotes 
an average over the jpd in (1). It is then natu- 
ral to consider the Laplace transform (e'^'P^'^^^) — 
JViG,N)e-'^P^^/'^dG, where the Laplace variable p is 
scaled by the factor f3N/2 for later convenience. In terms 
of the eigenvalues T^, this Laplace transform can then be 
expressed as a ratio of two partition functions 
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Using (1) 

ZpiN)=f ...[ exp[-f3E{{T,})]Y[dT, (4) 
Jo Jo 



where E {{T,}) = J^. {^T, - ^HT,)'^ 

In — Ti^\. Thus T^'s can be interpreted as 
the positions of charged particles confined in a 1-d box 
[0, 1], repelling each other via the 2-d Coulomb potential 
(logarithmic) and each subject to an external potential 
with a linear (with amplitude pN/2) and a logarithmic 
part, the latter being subdominant in the large iV ^ 1 
limit. Then E is the energy of a configuration and Zn{p) 
is the partition function at an inverse temperature p. 

Our next step is to evaluate the partition function 
Zp{N) in the large N limit. For that one carries out 
the multiple integral in (4) in two steps. The first step 
is a coarse-graining procedure where one sums over (par- 
tial tracing) all microscopic configurations of T^'s com- 
patible with a fixed charge density function gp{T) = 
J2i S{T — Ti) and the second step consists in per- 
forming a functional integral over all possible positive 
charge densities Qp{T) that are normalized to unity. Fi- 
nally the functional integral is carried out in the large N 
limit by the saddle point method. This procedure has re- 
cently been used in the context of the largest eigenvalue 
distribution of Gaussian [9] and Wishart random matri- 
ces [10] and also in other related problems of counting 
of stationary points in random Gaussian landscapes [11]. 
Following the general procedure in [9] , the resulting func- 



tional integral, to leading order in large N , becomes 

Zp{N) oc j I?[ep]e-^^'^['^-] (5) 
where the action is given by 

S[gp]^P f gp{T)TdT + B\[ gp{T)dT - 1 
Jo Uo 

- f f dTdT' gp{T)gp{T')\n\T -T'\. (6) 
Jo Jo 

Here _B is a Lagrange multiplier enforcing the normaliza- 
tion of gp{T). In the large N limit, one then evaluates 
the functional integral in (5) by the saddle point method, 
i.e., one finds the solution gp{T) (the equilibrium charge 
density that minimizes the action or the free energy) from 
the stationarity condition SS[gp]/dgp = which leads to 
an integral equation S = 2 g*{T') In \T-T'\ dT' . 
Differentiating once with respect to T leads to a singular 
integral equation 



(7) 



where Pr denotes the principal part. Assuming one can 
solve (7) for g^, one next evaluates the action S[gp] in 
(6) at the stationary solution and then takes the ratio 
in (3) to get 
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,-iN^-[S[g;]-S[gl]\^ (8) 



Inverting the Laplace transform gives the main asymp- 
totic result r{G,N) w cxp (^-fiV^^^ (-G^^ where 
the rate function is a Legendre transform, '^q{x) = 
maxp[— a;p -I- Jg{p)] with Jg{p) given by the free energy 
difference as in (8). 

It then rests to solve (7) to find the equilibrium charge 
density g^{T). Fortunately, singular integral equations 
such as (7) can be solved in closed form using Tricomi's 
theorem [12]. Skipping details [13], we find that the solu- 
tion g*{T) has three different forms (see Fig. 1) depend- 
ing on p, the parameter that controls the linear external 
potential. We get 
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Thus the equihbrium charge density undergoes two phase 
transitions as one tunes p respectively at p = 4 and p — 
—4. For p > 4, the Hnear potential is strong enough 
to confine the charges near T = leading to an inverse 
square root divergence of the density at T = and the 
density becomes zero beyond A/p. As p approaches the 
critical value 4 from above, the upper edge of the density 
support approaches the maximum value 1 and for —4 < 
p < 4 the charges are spread over the full box [0, 1] with 
the density diverging at the two end points. Finally for 
p < — 4, charges accumulate on the right wall T = 1 of 
the box with the density vanishing for T < 1 — 4/|p|. 
Note that the equilibrium density (in the large N limit) 
is independent of /?. The analytical solutions in the 3 
regimes above are depicted in Fig. 1 together with MC 
results [14]. 
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FIG. 1: Charge density Qp{T) for conductance (numerical, 
in red dots; theory, in solid black line). From left to right 
p = 6, —1, —6. The numerical density [11] is for N = 5 and 
(3 = 2. 

Inserting g*{T) in the action (6) one obtains the free 
energy difference in (8) 



Jg{p) 



'3/2+p + ln(-p/4) 

_Ei J- E 

32 2 

3/2 + ln(p/4) 



p < -4 

-4 < p < 4 (9) 
p > 4 



This analytical prediction is compared in Fig. 2 with MC 
simulations performed using similar tricks as in the case 
oi g*{T). Finally, the rate function ^'0(2^) = maxp[— a;p+ 
Jg(p)] in (2) is then given by (see Fig. 2): 
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Using the quadratic form in 1/4 < a; = G/N < 3/4 above 



in7'(G,iV) w cxp (^-f iV^^fG (^^jj = e "-v- 2; gives 

the Gaussian form in this central region from which one 
easily reads off the well known values for mean and vari- 
ance (G) = N/2 and var(G) = 1/8/3. On the other 
hand, near the two endpoints x = and x — 1, us- 
ing 'ifcix) ~ — ln(a;) and ~ — ln(l — x) as in (10), we 
get the power law dependence, V{G,N) - G'^^'/^ (^s 
G -> 0) and r{G,N) - (iV - G)''^'/^ (j^g c N) 
which are in agreement, to leading order in large N, with 



the exact far tails obtained in [S]. The central Gaus- 
sian regime matches smoothly with the two side regimes 
(only the third derivative is discontinuous). Such weak 
nonanaliticities of third order were also recently found 
in the entropy distribution in a quantum entanglement 
problem [15]. An expression similar to (10) was recently 
found in [7] by a different method. However, the inter- 
mediate regime with an exponential decay claimed in [7] 
does not appear in our solution. 

The same formalism can be easily extended to find the 
distribution V{P, N) of the shot noise P = J^^^i Tiil - 
Ti). It turns out to be convenient to work with the shifted 
variable Q = iV/4 - P = Ej=i(l/2 - Tif ■ Consider- 
ing the Laplace transform of the distribution of Q, we 
then get a new Coulomb gas where the charges, i.e., the 
shifted eigenvalues fjii = 1/2 — TiS are confined in the box 
[—1/2, 1/2], repel each other logarithmically and each ex- 
periences an external harmonic (with amplitude propor- 
tional to the Laplace variable p) -|-logarithmic potential. 
The equilibrium charge density Pp(/i) of this Coulomb 
gas can be computed again by solving the integral equa- 
tion p/i = Pr X^(/2 TT^^'^'^'' Skipping details [13], we 
again find two phase transitions upon tuning the value 
of p. For p > 8, the density has a semi-circular form over 
/i e {-^j2/p, \fVv\- For -8 < p < 8, the density has 
square root divergences at the box ends /j, = ±1/2. For 
p < —8, the density breaks up into two disjoint sectors 
around the two endpoints and vanishes in the middle 
of the box (see Fig. 3). In these three regions p > 8, 
—8 < p < 8 and p < —8 the exact form of the density is 
given by 
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where Cp = a/1/4 - 2/[p[ and clearly [Cp| < 1/2 for p < 
—8. In the last case, gp(/Lt) = for —C,p < < Qp. 

Following steps similar to the conductance case, we 
then compute the free energy difference associated with 
the shifted shot noise Q = N/A - P 



Jq{p) = 




(11) 



from which the rate function for the unshifted shot noise 
P = - Ti) can be computed. Noting that < 



4 



Ti < 1 implies that the scaled shot noise < x = P/N < 
1/4, we get (see Fig. 2): 
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Once again, the rate function ^'p(a; = P/N) is weakly 
nonanalytic at the critical values x = P/N = 1/16,3/16 
where a central Gaussian regime connects the two non- 
Gaussian regimes on either sides. From the central 
regime, one again reads off the mean and the vari- 
ance of the shot noise for large N: (P) = iV/8 and 
var(P) = 1/64/?. At the two endpoints P ^ and 
P N/A, we get using (12), power law behavior for 
the distribution, V{P,N) - pPN^/i (as P ^ 0) and 
V{P,N) ~ (7V/4 - P)^'^V4 (as P N/A). The latter 
limit is in agreement with the far tail result of [S], while 
the precise asymptotics near P ^ was not known be- 
fore. 
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FIG. 2: Left: Free energy difference Jg{p) (red) and Jq{p) 
(blue) from Montecarlo (MC) simulations and theory (solid 
line). We have used -30 < p < 30, /3 = 2 and iV = 5 (for 
conductance) and A'^ = 4 (for shot noise) respectively. Right: 
Rate functions ^G(a;) (green) and 'i'p{x) (blue) (see (10) and 
(12)). The black dots indicate the two critical points on each 
curve. 




FIG. 3: Charge density Qpi^i) for the shifted shot noise Q. 
From left to right p = 12, 2, -10. 

The large- A'' behavior (2) is generic for the distribu- 
tion of any linear statistics on the transmission eigenval- 
ues (of the form where a{T) can be nonlinear 
as in the case of shot noise) and the large-iV formal- 
ism developed here may, in principle, be extended to 
compute the distribution of any arbitrary linear statis- 
tic a(r). In particular, as a byproduct one can de- 
rive the mean and variance of any linear statistics by 



computing the central quadratic behavior of the asso- 
ciated rate function \E'(x). In certain cases this pro- 
vides new explicit formulae for mean and variance that 
are not easily accessible by other general methods such 
as in [16]. As an example we consider briefly the case 
of integer moments Tn — X^il^i^i"' relevant for statis- 
tics of charge cumulants [8, 17, 18]. The rate function 
close to the maximum (Gaussian regime) is given by [l'^]: 



^T,M = (4A„)-i(a;-P„)2 



for < X < X' 



(+) 



where x^"^ are n-dependent edge points and the con- 
stants An and P„ are given by yl„ = (2?! — 1)F(7t, + 
l/2)r(n- l/2)/167rnr2(n) and P„ = -i^'^Cn)- 



This implies: V{Tn,N) « exp 

near the mean. From this Gaussian form one can read 
off the mean and the variance. For the mean we find 
a /3-independent result (T;) = P„7V = ^(^^). In the 
special case (3 = 2 this was recently proved in [18] by 
other methods. For the variance we find a new exact re- 
suit: var(r„) = ^ = i (2"-i)r(n+i/2)r(»-i/2) ^^.^^ 

reduces to the conductance result 1/8/3 for n = I and 
approaches to the universal constant l/2/37r as n — > oo. 
Our new result for var(7^) for A'^ ^ 1 and its universal 
asymptote cannot be obtained easily from the formula in 
[18], valid for a finite number of open channels. 

In summary, we have obtained exact analytical dis- 
tributions for the conductance and the shot noise for a 
mesoscopic cavity with two leads and N channels each 
in the large N limit. Our results reveal a rich thermo- 
dynamic behavior including two phase transitions in an 
associated Goulomb gas problem leading to extraordinar- 
ily weak nonanaliticities in the distributions. We believe 
that the Coulomb gas method (well suited for large N) 
used here is very general and will be useful in other prob- 
lems as well. 
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